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ABSTRACT 

Kjeldsen et al. (1995) have developed a new technique for measuring steUar osciUations 
and claimed a detection in the G subgiant 77 Boo. The technique involves monitoring 
temperature fluctuations in a star via their effect on the equivalent width of Balmer 
lines. In this paper we use synthetic stellar spectra to investigate the temperature 
dependence of the Balmer lines, Ca 11, Fe i, the Mg b feature and the G band. We 
present a list of target stars likely to show solar-like oscillations and estimate their 
expected amplitudes. We also show that centre-to-limb variations in Balmer-line profiles 
allow one to detect oscillation modes with £ < A, which accounts for the detection by 
Kjeldsen et al. of modes with degree ^ = 3 in integrated sunlight. 

Key words: stars: oscillations - Sun: oscillations - Sun: atmosphere - stars: individual: 
ri Boo (HR 5235) - stars: individual: a Cen A (HR 5459) - S Scuti. 



1 INTRODUCTION 

, Measuring stellar oscillations is extremely difficult. The five- 
■ minute oscillations in the Sun, while rich in their information 
content, have tiny amplitudes. Observing similar oscillations 
in other stars therefore poses a great challenge (see Brown 
& Gilliland 1994 for a recent review). Most attempts have 
sought to detect periodic Doppler shifts of spectral lines. 
A second method involves measuring the total stellar lu- 
minosity, which varies due to changes in temperature in- 
duced by acoustic waves in the stellar atmosphere. Kjeldsen 
et al. (1995; hereafter referred to as KBVF) have proposed a 
new method in which they measure temperature changes via 
their effect on the equivalent widths of the Balmer hydro- 
gen lines. They presented strong evidence for oscillations in 
the G subgiant rj Boo (see also Bedding & Kjeldsen 1995), 
with frequency splittings that were later found to agree 
with theoretical models (Christensen-Dalsgaard et al. 1995a; 
Christensen-Dalsgaard, Bedding & Kjeldsen 1995b). 

Here we examine the equivalent-width method in more 
detail. In Section ^ we use synthetic spectra to calculate the 
temperature sensitivity of the Balmer lines and other strong 
absorption features. In Section |^ we use these results to esti- 
mate oscillation amplitudes for a list of target stars. In Sec- 
tion ^ we discuss the centre-to-limb behaviour of absorption 
lines and show that observing Balmer-line equivalent widths 
allows one to detect non-radial oscillation modes with £ < 4. 



2 TEMPERATURE DEPENDENCE OF 
EQUIVALENT WIDTHS 

The amplitude of an oscillation when measured in equivalent 
width (W) is related to the temperature fluctuations via 

SW _ dlnW ST 

W ^ dlnT T ' ^ ' 

Lines that are highly sensitive to temperature can be used 
to measure oscillations, while temperature-insensitive lines 
are useful as references. Also, detection sensitivity can be 
improved by combining observations of different lines, but 
only if the temperature dependence is known for each one. 

To determine the function W{T), we shall assume that 
an oscillating star can be represented at all times by a model 
atmosphere in radiative equilibrium. A more sophisticated 
approach would require that the acoustic waves be treated 
explicitly (see Breger 1975 and Frandsen 1984). With this 
assumption, we may calculate W{T) by considering a se- 
quence of normal stellar atmospheres spanning a range of 
effective temperatures. This is best done using theoretical 
models, since databases of observed stellar spectra are inho- 
mogeneous with respect to metallicity and surface gravity, 
and have uncertainties in T^s for each star. Note that in 
this paper we are not trying to estimate equivalent widths 
exactly. We are interested in relative changes as a function 
of effective temperature (and, in Section ^ of centre-to-limb 
position) . 
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Figure 1. Synthetic spectra of the G band and the Mg b feature 
for a model atmosphere having T^g = 5800 K. Equivalent widths 
were estimated by measuring the total flux in the regions enclosed 
by dashed lines. 



We have used the theoretical model GRS88 (Fuhrmann 
et al. 1993) to calculate spectra of Ha, H/3, Ca ii K, the 
Ca II H + He blend and the 4383 A Fe I line. Line profiles 
for the Balmer transitions were precalculated by K. Butler 
and T. Schoning (private communication) using VCS theory 
(Vidal et al. 1970). We have also calculated spectra of the 
G band and the Mg b feature (see Figure ^ using the LTE 
code RAIll by M. Spite, as extended for molecular lines by 
Barbuy (1982) — see Barbuy (1981, 1982) and Cayrel et 
al. (1991) for descriptions of the atomic and molecular data 
bases. The model atmospheres used in this case are from 
Gustafsson et al. (unpublished). 

For all models we have adopted solar surface gravity 
and metallicity; calculations of the Ha line by Fuhrmann et 
al. (1993) indicate that changing these parameters has lit- 
tle effect on equivalent widths. The microturbulent velocity 
was assumed to be 1 km/s (note this value is not critical 
for the strong lines we are considering) and the damping 
constants were obtained by fitting to the solar spectra (and 
also the Arcturus spectrum, for the G band and Mg b). 
The models do not include the chromospheric temperature 
inversion, but chromospheric effects should have negliglible 
influence on the equivalent widths of the metal lines consid- 
ered here because only the innermost line cores are formed 
at logrsooo < —3.5. For the Balmer line, equivalent widths 
near the limb are signiflcantly influenced by non-LTE effects 
in the line core, as discussed in Section [O below. 

The upper panel in Figure |^ shows W{Tas) as measured 
from our synthetic spectra. The lower panel shows the slope 
{d\nW / d\nTc.s), which we obtained simply by measuring 
the differences between adjacent points in the upper panel. 
The slope for the Balmer lines is 5-7, consistent with the 
value of 6 adopted by KBVF, which was based on published 
models (Kurucz 1979; Gray 1992). We see that the Mg b 
feature is especially promising for oscillation measurements. 
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Figure 2. Equivalent width and slope (9 In In Teff ) for var- 
ious lines as a function of stellar effective temperature, based on 
model atmosphere calculations. 



The Fe I lines are also useful since, although they are in- 
dividually weaker, they are numerous and very sensitive to 
temperature. 



3 OSCILLATIONS IN TARGET STARS 

In the Sun, the amplitudes of individual oscillation modes 
can vary considerably over timescales of several days. How- 
ever, different modes vary independently and the maximum 
amplitude, taken over all modes, stays roughly constant. 
When the solar oscillations are measured in bolometric lu- 
minosity, this amplitude is ((5L/L)boi = 4.1 ppm (parts-per- 
million) — see Kjeldsen & Bedding (1995; hereafter referred 
to as KB) and references therein. 

From KB, equations 5 and 7, amplitudes of solar- like 
oscillations in other stars should scale roughly as 



(5L/L)boi = 



(M/M© )(reff/5777K) 



4.1 ppm. 



(2) 



This luminosity variation is due almost entirely to changes in 
temperature (the change in radius is negligible). Therefore, 
using L oc R^T"^ and Equation hi we can write 



5W dlnW 



1.0 ppm. 



(3) 



W ainTeff (M/M©) (reff/5777K) 

Equation |^ gives an estimate for the oscillation signal 
expected from equivalent-width measurements of a spectral 
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Table 1. Expected oscillations signal from observations of Ho equivalent width for a sample of bright stars 
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density, as described below. We also omitted known vari- 
ables (e.g., S Scuti stars) and close binaries having nearly- 
equal components. For the remaining binaries we give the 
parameters for the primary only (except a Cen, whose com- 
ponents are listed separately). 

For the Sun, a Cen A and B, Procyon (a CMi), rj Boo, 
f3 Hyi and e Eri we adopted the same fundamental param- 
eters as used in KB and KBVF. For the other stars, we 
estimated effective temperatures using T^s — 11000 K/(_B — 
V + 1.24), used parallaxes from The Bright Star Catalogue 
and took bolometric corrections from Lang (1992). We made 
no reddening corrections except for HR 1543, which has 
E{B — V)— 0.03. Masses were estimated using stellar mod- 
els in Fig. 1 of Bressan et al. (1993), which assumes solar 
metallicity, OPAL opacities and convective overshoot. 



line. The photon noise, on the other hand, scales as the 
square root of the number of photons detected in the line. 
Hence, assuming photon noise to be the dominant noise 
source, the signal-to-noise ratio for a fixed observing time 
will scale as: 



S N (X — — , 4 

d In Toff MTcB 

where m is the stellar magnitude at the relevant wavelength. 

We have applied Equations |^ and 4 to a sample of bright 
stars that are likely to undergo solar-like oscillations. Ta- 
ble ^ shows fundamental parameters and expected oscilla- 
tion characteristics for the stars, which were selected from 
The Bright Star Catalogue (Hoffleit 1982) to have V < 4.2 
and 0.35 < _B — V < 1.2. A further restriction to include only 
main-sequence stars and subgiants was made on the basis of 
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For each star we used Equation ^ to calculate the oscil- 
lation signal expected from observations of equivalent- width 
fluctuations of Hq. Computing similar results for other lines 
is straightforward using the data in Figure ^. The predic- 
tions for 1] Boo and the Sun agree roughly with the Balmer- 
line observations of KBVF, although the amplitude scale in 
Table |^ may need revision once more observations become 
available (we are currently analysing high-quality data from 
Q Cen A, Procyon and the Sun). However, the relative am- 
phtudes of the stars should not change drastically unless 
Equation |^ which is based on theoretical models, is found 
to be incorrect. 

The stars in Table ^ are ordered according to S/N. This 
is the signal-to-noise ratio for a fixed observing time, arbi- 
trarily normalized to give S/N — 3 for the Sun at a distance 
of 1 pc. How much observing time would this require? Based 
on the KBVS observations of 77 Boo, a total of 1.4 x 10^° 
photons/ A are required to produce a photon-noise level of 
lOppm. This allows us to estimate the observing time re- 
quired to achieve a given S/N. For example, a 2.5-metre 
telescope collecting photons with a total efficiency of 5% 
would take about 43 hours to achieve the S/N given for each 
star in Table |l[ This estimate is based on observations at 
Ha only - simultaneous measurement of other spectral lines 
would improve matters. Also, we have neglected 1/f noise 
sources, such as instrumental drift, which will be important 
for stars oscillating at low frequencies (i.e., sub-giants). 

The ranking of targets in Table ^ may differ for the other 
observing techniques (luminosity and velocity) and can be 
calculated using the formulae in KB. The Table also includes 
rough estimates for the frequency of maximum oscillation 
amplitude fmax and the large frequency separation A^o, 
computed by scaling from the Sun (see KB) . These frequen- 
cies will be important when choosing targets and designing 
observing programs. The large separation depends on the 
average density of the star and the table only includes stars 
for which Avo > l/(day) = 11.57/^Hz, in order to restrict 
the sample to main-sequence stars and subgiants. 

The table also shows the quantity 



Urot sm I 



V smi/{2nR) 

V sini/kms"'^ 



0.23 ^Hz- 



(5) 
(6) 



where Vrot is rotation frequency of the star. Components of 
an oscillation mode that is split by stellar rotation will be 
separated by approximately i/rot (see below). 

Finally, it is worth mentioning the double-lined spectro- 
scopic binary a Equ (HR 8131), whose components (G5 III 
and A5 V) are slightly too faint to be included in Table |l| 
The masses and luminosities of both components have been 
derived from interferometric and spectroscopic observations 
by Armstrong et al. (1992), making this system particularly 
valuable for testing stellar models. 



4 MODE SENSITIVITY 
4.1 Calculations 

From equivalent-width observations of the daytime sky, 
KBVF detected solar oscillations with degrees ^ = 0, 1, 2 
and 3. The detection of / = 3 was puzzling, given that the 
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Figure 3. Observed Ha spectra of the Sun at the centre and edge 
of the disk, taken from Brault &; Testerman (1972). The wings of 
the Ha line are much weaker in the lower diagram, whereas the 
metal lines and the core of Ha show little variation. 



the observations did not resolve the solar disk. The answer 
is straightforward: the Balmer lines are much weaker at the 
edge of the solar disk than at the centre (see Figure ^), so 
their integrated profiles are strongly weighted towards the 
disk centre. This has lead us to investigate in detail the 
sensitivity of the equivalent-width method to modes with 
different values of £. 

The calculation involves integrating the oscillations over 
the stellar surface, including the effects of limb darkening 
(Christensen-Dalsgaard & Cough 1982; Heynderickx et al. 
1994; Robinson et al. 1995). The oscillations of a star can be 
described in terms of spherical harmonics Y[^{6,(j)). If the 
star is rotating slowly, there is no preferred direction and we 
are free to choose a coordinate system with the polar axis 
pointed towards the observer. In this case, all integrations 
of spherical harmonics will be zero unless m = Q. Even for 
stars in which rotation is important, the following analysis 
is still useful, as discussed at the end of this section. 

For m = 0, the spatial response function for each I (the 
ratio of the observed to the actual amplitude) is 



Se = 2V2l + l Pe{fi)Iifj,)F{fi)fidf^. 



(7) 



Here /j. = cos 61, Pei^y) is the Legendre polynomial of de- 
gree I, /(/i) is the centre-to-limb variation of continuum in- 
tensity (classical limb darkening) and F{fi) is the centre-to- 
limb variation of the oscillation signal. 

The function /(/i) for the Sun is well approximated by 



/(At) = l-tt2(l-M)-t'2(l 



M ), 



(8) 



where the coefficients vary with wavelength (Allen 1973). 
The function F{fj,) depends on how the observations are 
made. If the oscillations are observed via the Doppler shift 
of a spectral line then F{^) = ^t, corresponding to the pro- 
jection of a nearly-radial pulsation along the line of sight. 
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Figure 4. Centre-to-limb variation of equivalent width from 
solar-temperature model atmospheres, normalized to unity at the 
centre of the disk. 



For observations in intensity there is no projection effect (to 
a good approximation) and F{iJ,) — 1. If the oscillations are 
observed in equivalent width then we set F{fi) — the 
spatial variation of equivalent width from centre to limb. 

In principle, we could determine the function W{fi) 
from spatially-resolved spectra of the Sun. However, at visi- 
ble wavelengths it seems that only two have been published: 
the Preliminary Kitt Peak Photoelectric Atlas (Brault & 
Testerman 1972), which has spectra at /i = 1.0 and /i = 0.2 
(see Figure |^, and the atlas by Delbouille et al. (1973), 
which contains a single spectrum at = 1.0. 

To estimate W{^) at all values of /x, we have therefore 
used the models described in Section ^ to produce intensity 
spectra at different points on a stellar disk. The results for a 
star with solar temperature, normalized to the disk centre, 
are shown in Figure ^ For the Balmer line cores, non-LTE 
effects are important and the models significantly underesti- 
mated the area of the line cores near the limb. To account for 
this, we assumed the difference between the actual equiva- 
lent width and that predicted by the GRS88 model increases 
linearly with decreasing /i. Based on the measurements of 
VK(/i) at two points on the solar disk (/^ = 1.0 and 0.2), we 
could correct the GRS88 equivalent widths at all intermedi- 
ate points (the correction at = 0.3 was 35% for Ha and 
18% for H/3). The final values are shown in Figure ^ 

For those lines that were computed using the GRS88 
model, we have repeated the calculations using the semi- 
empirical model of Holweger & MuUer (1974). For the 
Balmer lines, we again corrected the equivalent widths for 
the missing core area. The resulting W(p) profiles for all 
lines are similar to those in Figure 111 to within 10%. 

One other comparison with observation is possible. 
Quoting the unpublished Michigan Atlas, Holweger (1967) 
gives equivalent widths at ^ = 1.0 and /i = 0.3 for many 
lines, including 5172 A Mg I and 4383 A Fe I (but note that 
his values for Ca ll H & K and H I come from model calcu- 
lations). For these lines, the ratio W(/i — Q.3)/W(/i — 1.0) 
given by Holweger agrees with our calculations to within 



10%. It would clearly be desirable to have detailed centre- 
to-limb measurements of all these strong lines, similar to 
those of the 7774 A O I triplet by Affrock (1968) and King 
& Boesgaard (1995). 

Our general conclusion is that the Balmer lines are 
strongly limb darkened, while the other lines have almost 
constant strength from centre to limb. For our purposes, it 
is sufficient to approximate F{^) for each spectral line using 
a linear function: 



F{^l) = l-c{l-^J). 



(9) 



By fitting a straight line to each set of points in Figure ^ 
(giving more weight to the central portion of the disk, since 
this contributes the greatest area), we obtained the following 
values for c: 0.95 for Ha, 1.0 for H/9 (which we also adopt 
for H7 and H5), -0.2 for Ca 11 K, 0.1 for the Ca 11 H + He 
blend, -0.2 for Fe I, 0.05 for the G band and 0.15 for the 
Mg b feature. 

Our Balmer-line models indicate that c varies quite 
slowly with effective temperature (c oc T^'^), so the results 
given here apply to a broad class of stars. 

Having obtained approximations for 7(^) and F{ii), we 
can now integrate Equation ^. Results for the five lowest- 
degree modes {£ — 0, . . . , 4) are: 
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(10) 



This matrix product, which should be evaluated from right 
to left, can be used to estimate the mode sensitivity of each 
observing method as a function of wavelength. Note that 
Equation 10 can be used for velocity measurements by set- 
ting c = 1 and for intensity measurements by setting c = 0. 



4.2 Results and discussion 

Figure |^ shows the sensitivities as calculated from Equa- 
tion 10, expressed relative to ^ = 0. In the case of velocity 
and intensity observations (solid and dotted curves), there 
are published solar oscillation data with which to compare 
the calculations: 

(i) Observations in velocity from the South pole using the 
Na I line (590 nm) by Grec et al. (1980) give Si/So = 1.5, 
S2/S0 = 1.2, S3/S0 = 0.6 and S4/S'o = 0.2 (the scatter is 
about ±0.2). 

(ii) Observations in bolometric intensity from the So- 
lar Maximum Mission by Woodard & Hudson (1983) give 
Si/So = 1.25 and S2/S0 = 0.83 (scatter about ±0.06). 

(iii) Observations in luminosity with the IPHIR green 
channel (500 nm) by Toutain & Frohlich (1992) give Si/ So = 
1.09 and S2/S0 = 0.79 (scatter about ±0.07). 

In each case, there is good agreement with our results. 
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Figure 5. Sensitivity of observing methods as a function of wave- 
length. The solid line is for observations in velocity and the dashed 
line for intensity. Individual points indicate equivalent-width ob- 
servations in different spectral features, with symbols having the 
same meaning as in Figures ^ and ^ (but with the addition of If7 
and H(5, shown as open diamonds). Negative values of Si mean 
the oscillations will appear to have reversed phases. 



Turning to equivalent-width measurements, we see from 
Figure |5| that the Balmer lines give a response similar to 
that of velocity measurements. This follows from their strong 
centre-to-limb variation and is consistent with the solar ob- 
servations by KBVF — in particular, we can now explain 
their detection of ^ = 3. The other absorption lines consid- 
ered here have much weaker centre-to-limb variations (c ~ 0; 



see Figure y) and so have spatial responses similar to that 
of intensity measurements. This suggests the possibility of 
determining the I value of a given mode using simultane- 
ous observations in several absorption lines, which might be 
interesting for studies of 5 Scuti variables (Matthews 1993). 

Finally, we consider a rotating star. A mode with a 
particular (n,£) will split into a multiplet having m = 
— 1,1, with the components being separated by 
about v-cot (the oscillation power summed over these peaks 
is conserved). If the star is seen pole-on, only modes with 
m = are observable and the analysis of this section holds 
exactly: we would not notice that the star is rotating. At 
the other extreme, in a star seen from near the equator 
(e.g., the Sun) the modes with \l — m\ odd will have zero 
amplitude because the Legendre function is antisymmetric 
around the equator. For ^ = 1 we would therefore only see 
two modes (m = ±1). The total power will be the same as 
for a non-rotating star, so the amplitudes of these two modes 
will each be multiplied by l/\/2 relative to the non-rotating 
case. For £ = 2 we would see three modes (m = —2, and 2) 
with relative amplitudes of -^^3/2, 1/2 and yJ'i/2. Other 
cases can readily be calculated using the formulae given by 
Christensen-Dalsgaard (1994). 
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